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Abstract. - Quasi-one-dimensional superconductors or nanowires exhibit a transition into a 
nonsuperconducting regime, as their diameter shrinks. We present measurements on ultrashort 
nanowires (~40-190 nm long) in the vicinity of this quantum transition. Properties of all wires 
in the superconducting phase, even those close to the transition, can be explained in terms of 
thermally activated phase slips. The behavior of nanowires in the nonsuperconducting phase 
agrees with the theories of the Coulomb blockade of coherent transport through mesoscopic 
normal metal conductors. Thus it is concluded that the quantum transition occurs between 
two phases: a "true superconducting phase" and an "insulating phase". No intermediate, 
"metallic" phase was found. 



Under certain conditions, usually associated with a critical resistance per square [1,2], crit- 
ical total resistance [3-6], or a characteristic diameter [7,8], a wire made of a superconducting 
metal looses its superconductivity and acquires two signatures of insulating behavior: i) The 
resistance increasing with cooling and ii) a zero-bias resistance peak [3,4]. There are many 
models that capture certain features of the SIT in ID wires. Some rely on the "fermionic" 
mechanism, in which disorder combined with electron-electron repulsion suppresses the crit- 
ical temperature, Tc, to zero [9]. In other, "bosonic", models the order parameter remains 
nonzero in the "insulating" phase while the coherence is destroyed by proliferating quan- 
tum phase slips (QPS) [10-15]. Existing theoretical models frequently predict a quantum 
superconductor-insulator transition (SIT) in thin wires [11,13,16,17], driven, in many cases, 
by the interaction of the fluctuating phase with the Caldeira-Leggett environment [18]. Con- 
ditions that make QPS experimentally observable and the relation of the QPS to the SIT are 
still being actively researched [1-8, 19-22]. 

Here we present a quantitative analysis of the transport properties of ultrashort nanowires 
in each of the two phases - the insulating phase and the superconducting phase. We show 
that the insulating phase is characterized by the normal-electron transport and governed by 
the Coulomb blockade physics [23,24]. The wires in the superconducting phase exhibit good 
agreement with the Langer-Ambegaokar-McCumber-Halperin (LAMH) theory of thermally 
activated phase slips (TAPS) [25-27], without any QPS contribution. The TAPS physics is 
dominant, even in the vicinity of the SIT. Thus we conclude that the observed transition 
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Fig. 1 - (a) SEM micrograph of an ~8 nm wide nanowire (light) suspended over the trench (black) in 
SiN. The sputtered MoGe film was 5.5 nm thick. The white regions at the ends of the wire indicate 
that this wire is suspended straight, without kinking and entering into the trench. The R{T) curves 
for insulating and superconducting wires are shown in (b) and (c). The arrow in (b) shows Rn for 
sample D. In (c), solid curves indicate fits to the LAMH-TAPS theory. The fitting parameters are 
the coherence lengths, 70.0, 19.0, 11.5, 9.4, 5.6, and 6.7 nm, and the critical temperatures, 1.72, 2.28, 
3.75, 3.86, 3.80, and 4.80 K, for samples 1-6, respectively. The corresponding normal resistances and 
the lengths, determined from the SEM images, are 5.46, 3.62, 2.78, 3.59, 4.29, 2.39 kfi and 177, 43, 
63, 93, 187, 99 nm respectively. 



occurs between a truly superconducting phase (which shows no QPS and thus the resistance 
approaches zero resistance at T = 0) and an insulating phase in which the wire is in the 
normal state and the transport is controlled by weak Coulomb blockade. 

The nanowires, ranging in length between 43 and 187 nm, are fabricated by sputtering 
of amorphous Mo79Ge2i alloy on top of suspended fluorinated single-wall carbon nanotubes 
[3,4,20]. The wires are homogeneous as is seen from scanning electron microscope (SEM) 
images (Fig. [T^). The electrodes are deposited during the same sputtering run as the wire 
itself. Since Ar-atmosphere sputter-deposition is isotropic, and due to the small diameter of 
the nanotubes (~l-2 nm), the wires, which occur on the outer surface of the nanotube, form 
seamless connections to the electrodes. The homogeneity of wires is confirmed also by the fact 
that their normal resistance is close to that estimated from the sample geometry and known 
bulk resistivity (~200 /j,r2-cm) [3,28]. Transport measurements are performed in "^He and "^He 
cryostats equipped with leads filtered against electromagnetic noise [29] . One sample (sample 
F) was measured down to ^^20 mK. 

Resistance vs. temperature, R(T), data for insulating and superconducting samples are 
shown in Figs. [T]d and [TJ;, respectively, where R{T) = dV/dl a,t V 0. The resistive 
transition observed in all samples at higher temperatures is that of the film electrodes, which 
are connected in series with the wire. As the electrodes go superconducting, the total sample 
resistance equals the wire resistance. The only changing parameter amongst all samples is 
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the nominal thickness of MoGe (4.0-8.5 nm). Consistently, the critical temperature of the 
electrodes decreases gradually as we proceed from the superconducting sample corresponding 
to the largest amount of MoGe sputtered (sample 6) to the insulating sample with the smallest 
amount of MoGe sputtered (sample B) . The resistance measured immediately below the film 
transition is taken as the normal (or high-temperature) resistance of the wire, Rn- 

The R{T) curves of superconducting samples are shown in Fig.[TJ:. The fits are made using 
the LAMH-TAPS formulas [20]. As the resistance starts to sharply drop with the cooling, the 
LAMH model becomes valid, and the data exhibit an excellent agreement with the fits. For the 
two thinnest samples (1 and 2), which have the lowest T^'s, the thermodynamic critical field 
Hc{T) in the free energy barrier for phase slips had to be modified: We used the empirical 
expression [30] Hc{T) oc 1 - [T/T^f, instead of the usual Hc{T) oc 1 - T/T^. After such 
modification the fits matched the data. A striking result is that our set of ultrashort wires 
exhibits an excellent agreement with the TAPS model, without requiring any QPS contribution 
[7, 10]. This is even true for wires in the vicinity of the SIT. Good agreement with the LAMH 
model indicates that the observed superconducting regime is a "true" superconducting phase, 
i.e. the wire resistance is expected to rapidly approach zero as T — > 0. In this regime, the 
wires R{T) exhibits a negative curvature on ln(i?) vs. T plots, which is an indication that the 
contribution of QPS [7, 10, 19,21,31] is negligible. Such QPS-free regime is new and was not 
seen on longer wires [7,8, 19,22]. 

The R{T) curves of nonsuperconducting wires show a qualitatively different behavior 
(Fig- [lb)- We term them "insulating" because they reproducibly show (i) an upturn at the 
lowest temperatures, i.e. dR{T)/dT < (down to ^20 mK, as was tested for sample F), and 
(ii) a zero-bias resistance peak, i.e. (PV{I)/dP < for ^ (/ is the bias current and V is 
the bias voltage). 

The observed abrupt change from TAPS to the insulating behavior strongly suggests that 
a quantum phase transition does occur in ultrashort nanowires. This is in contrast to longer 
wires, which exhibit a crossover [7,8,19,22] from a quasi-superconducting to a quasi-normal 
regime. We speculate that this SIT takes place due to coupling of QPS to gapless excitations 
in the environment [11,13-16], similar to the Chakravarty-Schmid transition [31-34]. The 
QPS-free regime can be understood assuming that the QPS are completely suppressed by a 
Caldeira-Leggett environment (e.g. produced collectively by the QPS cores). In the insulating 
phase the QPS proliferate and completely suppress superconductivity, again due to normal 
cores associated with each QPS. 

To understand dR{T)/dT < seen in our insulating samples, we consider the theories of 
the Coulomb blockade (CB) in diffusive normal wires [23,24,35,36]. Nazarov showed that 
the CB can survive in a setting in which two plates of a capacitor C are connected by a 
homogeneous normal wire (which now plays the role of a tunnel barrier), even if its resistance 
is much lower than the von Klitzing constant, Rk — h/e'^, provided that the wire acts as a 
coherent scatter [23]. Golubev and Zaikin (GZ) [24] derived useful I{V) formulas, enabling 
direct comparisons with experiments. At high temperatures [ksT > Ec , where Eq = e^/2C 
is the charging energy) the zero-bias conductance, G(T) = 1/R{T) is 



where (3 = 1/3 for diffusive wires. Also, C(3) = 1.202, g = GqRk, and Go is the conductance 
in the absence of the CB. Apart from the second order term {Ec/ksT)'^ , this result is the 
same as obtained by Kaupinnen and Pekola for a single tunnel junction [37]. Originally such 
result was derived for a primary thermometer [38] . The same expression was derived by Joyez 
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Fig. 2 - Pekola et al. [38] conductance function 7(T) = Go [Go — G(T)] ^ plotted vs. temperature 
for all insulating samples. The solid lines are fits to the GZ theory (Eq. (O). Values of the fitting 
parameter Go for samples A-H are Gq ^ =6.14, 7.14, 7.93, 7.76, 9.78, 9.97, 17.33, and 26.10 kQ,, 
respectively. As expected, they are close to the corresponding normal resistances of the wires: _Rjv = 
6.43, 7.54, 8.25, 8.35, 10.33, 10.50, 18.05, and 32.46 kVL. The predicted range of the offsets is shown 
by arrows. The lengths of these wires are 46, 45, 140, 105, 140, 49, 120, and 86 nm. 

and Esteve [39] for a single tunnel junction (i.e. for the dynamic Coulomb blockade), though 
in their case the value of parameter g is defined differently: Gje = RK/Renv {Rem is the 
impedance of the environment). If the ratio Ec/ksT is a small parameter, the diverging 
terms of the Eq. ^ can be removed by rewriting it as: 



In this form the approximate Eq. ^ follows closely (as our numerical analysis shows) the 
exact result of GZ given in Eq. (28) in Rcf. [24]. In Fig. [2]we rcplot the R{T) data according 
to Eq. ([2]). The predicted linearity of the Pekola function 7(r) is indeed observed on all 
insulating samples, confirming that CB does occur. The parameter Gq, adjusted to produce 
the best linearity of the curves, matches the high-bias conductance, as expected (see the 
caption to Fig. The Ec (Eq. ^) is determined from the slopes. As indicated by Eq. 
the offset at T = should be larger than zero. The arrows in Fig. [5] show the expected 
range of the offsets, with g = Rk/Rn, which is in a fair agreement with the experiment. At 
higher temperatures (for example at T > 1 K for sample A) we observe a deviation from the 
linear dependence, which might be due to weakening of the proximity effect, induced by the 
superconducting leads. 

At higher bias, the CB appears as a zero-bias anomaly (ZBA) observed in the differential 
conductance versus bias voltage dependence for all our insulating wires. The anomaly is 
similar to the one measured by Pekola et al. on single-electron transistors [38]. Examples of 
ZBA are shown in Fig. [3] We compare the shape of the dI{V)/dV curves to the GZ theory, 
in which the ZBA in the I{V) dependence is given as 



where w = u -\- iv, u = gEc/T^^ksT, v — eV/TrksT, and '^{x) is the digamma function. The 
fitting was done by differentiating Eq. ^ and making an expansion of '^{x). The results 
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Fig. 3 - (a) The dl/dV vs. V curves for sample B at T = 0.28 K . Comparisons to the GZ theory 
are shown with q^ff — e (dashed hne), q^ff = 1.29e (sohd hne), and qeff ~ 2e (dotted hne). (b) The 
dl/dV vs. V curves for sample E at T = 0.3 (deepest dip), 0.5, 0.75, 1.0, and 1.5 K (shallowest dip). 
q^ff — 1.53e. (c) The dl/dV vs. V curves for sample D at two different magnetic fields. At B = 
the leads are superconducting and qeff = 1.70e whereas at high field (B = 9 T) the leads are driven 
normal and q^ff drops to 1.05e. Solid lines are fits to the GZ theory. 



obtained on high i?jv samples agree well with this prediction. On the other hand, for samples 
near the SIT point (i.e. Rn ~ h/Ae^ [31]), the width of the conductance dip was narrower 
than predicted (Fig. [3^). Since theoretically the resistance peak becomes narrower as f ~ e/T 
increases, the narrowing of the ZBA cannot be explained by an electronic heating, as this would 
lead to a wider anomaly. Understanding this peak-narrowing effect is difficult because the 
theory is derived for normal electrodes. We find empirically then that the dI{V)/dV obtained 
from Eq. ^ can be used to fit our data with e replaced by an effective charge, Qeff, which we 
used as a fitting parameter along with Go and Ec- With this correction the model matches 
the experiments. We speculate that the Andreev reflection taking place at the ends of the wire 
leads to Qeff > e. We also verified that the dliV) / dV curves measured at various temperatures 
show a good agreement with theory (Fig.lSh)- Finally, if a perpendicular magnetic field is used 
to drive the electrodes normal (Fig.[3J;), the effective charge qeff decreases to e, as expected for 
all- normal systems. Another explanation for the peak narrowing is the refrigeration effect [40]. 

In Figs. and [IJa we compare the parameters Gq and Ec obtained from the fits to 
the dI{V)/dV plots (measured at 0.3 K) and the same parameters extracted from the fits 
to the 7(T) data. A close agreement is observed, thus confirming the consistency of our 
interpretation. Furthermore, it is clear from Fig. [4^ that Go is very close to R^^ for all 
samples, as it should be. 

The Ec extracted from fitting to the ^{T) and the dI{V)/dV plots is shown in Fig. |3|d as 
a function of R^. A strong decrease of Ec with decreasing Rjy is seen, possibly due to the 
renormalization with Rn [23] . The exact physical meaning of the Ec is not well established. 
One possibility is that the dynamic Coulomb blockade [39] that involves just one barrier (i.e. 
the entire nanowire) is realized and the effective capacitance is the capacitance between the 
electrodes (i.e. two coplanar thin films [3]). The relevant sections of the electrodes, which 
define the effective value of G, are determined by a "horizon" [37,41]. To roughly estimate 
the charging energy (Ec = e^/2G) we calculate the capacitance of the two coplanar plates 
as G = eeQwK{k)/K{l — fc^)^/^, where e « 2 and K{k) is the complete elliptic integral of 
modulus fc = (1 + L/a)^^ [42]. For our electrodes of width w = 10 fim, the length a = 50 
/im , and the spacing between closest edges L = 0.1 /xm, the capacitance is estimated G « 1 
fF and the corresponding charging energy is E^max ~ 80 iieV (top dashed hne in Fig. HJ^). 
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Fig. 4 - (a) Values of Go, extracted from the GZ theory fits, plotted versus ■ Symbols are 
Go obtained from temperature-dependent (x) and voltage-dependent (o) data. The dashed line 
represents Go = R]^^ ■ (b) Charging energies extracted from fits to the GZ theory and plotted versus 
Rn- Symbols are Ec obtained from the temperature-dependent (x) and voltage-dependent (o) data. 
Dashed lines indicate the range of Ec estimated for the geometry of our electrodes. 



The lower limit, Ec.min, was obtained taking into account the capacitance corresponding to 
the entire area of the electrodes, including the contact pads. Using the program FASTCAP 
(www.fastfieldsolvers.com) we obtained Ec,min ~ 8 eV (lower dashed line in Fig. Most 
of the experimental Ec values fall between the two limits (Fig.^jD). However, since the length 
of the studied wires is roughly of the same order as the dephasing length L^, we note that 
for wires with L > there is another model [43], which predicts the same form of 7(T) 
as Eq. ([2]) but uses an effective charging energy that depends only on the wire parameters: 
Ec = 3/27r27Vod^ [44], where No ~ 8 x lO^^ (eV cm)-^ is the density of states of MoGe and 
d is the diameter of a wire. This model gives roughly the same values of the charging energy. 
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